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I. INTRODUCTION 



Meson-antimeson mixing is an important vehicle for indirect studies of New Physics (NP). 
Due to the absence of tree-level flavor- changing neutral current transitions in the Standard 
Model (SM), it can only occur via quantum effects associated with the SM and NP particles. 
In fact, the existence of both charm and top quark were inferred from the kaon and beauty 



mixing amplitudes The estimates of masses of those particles were later found to be in 
agreement with direct observations. This motivates indirect searches for NP particles in a 
meson-antimeson mixing. 

Recently, there has been a considerable interest in the only available meson-antimeson 
mixing in the up-quark sector, the D° — D° mixing . The fact that the search is indirect 
and complimentary to existing constraints from the bottom-quark sector actually provides 
parameter space constraints for a large variety of NP models [3j, |4|. 

A flurry of recent experimental activity in that field led to the observation of D° — D° 
mixing from several different experiments such as BaBar fl, Belle and CDF Q. These 
results have been combined by the Flavor Averaging Group (HFAG) [8] to yield 

y e D xp = (6.6 ±2.1)- 1(T 3 (1.1) 

x%p = 8.7± 3 ;° • 10" 3 , (1.2) 

where xd and yn are defined as 

- AMd a AFd h *\ 

xd = —= — , and y D = — — , (1.3) 

Id 

Td is the average width of the two neutral D meson mass eigenstates, and AMd, AFd are 
the mass and width differences of the neutral D-meson mass eigenstates. In the limit of 
CP- conservation, Ar^ = T + — r_, where "+" and "-" are CP-even and CP-odd D-meson 
eigenstates respectively. 

One can also write yr> as an absorptive part of the D° — D° mixing matrix 9|, 



r D 



n 



where p n is a phase space function that corresponds to a charmless intermediate state n. 
This relation shows that Ar^j is driven by transitions D°,D° — > n, i.e. physics of the 
AC = 1 sector. 
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Eqs. (11. ip and (11.21) imply one-sigma window for the HFAG values of xd and y-o, 

5.3- 1(T 3 <x D < 11.7- 1(T 3 (one — sigma window) (1-5) 
4.5- 1(T 3 < y D < 8.7- 1(T 3 (one — sigma window) (1-6) 

In principle, these results can be used to constrain parameters of NP models with the 
anticipated improved accuracy for the future D-mixing measurements. In reality, those 
results can only provide the ballpark estimate to be used for constraining NP models. The 
reason is that the SM estimate for the parameters xd and y^> is rather uncertain, as it is 
dominated by long-distance QCD effects It was nevertheless shown that even this 

estimate provides rather stringent constraints on the NP parameter space for many models 
affecting the mass difference x& {3}, jl^-fisl]. 



It was recently shown {4] that D° — D° mixing is a rather unique system, where y^> can 
also be used to constrain the models of New Physics 1 . This stems from the fact that there is 
a well-defined theoretical limit (the flavor £[/(3)-limit) where the SM contribution vanishes 
and the lifetime difference is dominated by the NP AC = 1 contributions. In real world, 
flavor SU(3) is, of course, broken, so the SM contribution is proportional to a (second) 
power of m s /A, which is a rather small number. If the NP contribution to yv is non-zero in 
the flavor S77(3)-limit, it can provide a large contribution to the mixing amplitude. 

To see this, consider a D° decay amplitude which includes a small NP contribution, 
A[D° -> n] = + A^ p \ Experimental data for D -meson decays are known to be in a 



decent agreement with the SM estimates 20, |2]J . Thus, A^ p ^ should be smaller than (in 



sum) the current theoretical and experimental uncertainties in predictions for these decays. 
One may rewrite equation (II. Ah in the form (neglecting the effects of CP-violation) 

yn = Y, ^4 SM) 4 SM) + 2 E ^4 NP) 4 SM) + E ^4 NP) 4 NP) • (1.7) 

n 1 D „ J- D n 1 D 

The first term in this equation corresponds to the SM contribution, which vanishes in the 
SU(3) limit. In ref. [4j the last term in (11.71) has been neglected, thus the NP contribution 
to 2/d comes there solely from the second term, due to interference of A^ M ^ and A^ p \ While 
this contribution is in general non-zero in the flavor SU (3) limit, in a large class of (popular) 



models it actually is 



221 ]. Then, in this limit, yn is completely dominated by pure A^ F ^ 



1 A similar effect is possible in the bottom-quark sector 19 1. 
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contribution given by the last term in eq. (11. 7p ! It is clear that the last term in equation 
(I1.7P needs more detailed and careful studies, at least within some of the NP models. 

Indeed, in reality, flavor 577(3) symmetry is broken, so the first term in Eq. (jl.7p is not 
zero. It has been argued [h]] that in fact the SM SU (3)-violating contributions could be at 
a percent level, dominating the experimental result. The SM predictions of yo, stemming 
from evaluations of long-distance hadronic contributions, are rather uncertain. While this 
precludes us from placing explicit constraints on parameters of NP models, it has been argued 
that, even in this situation, an upper bound on the NP contributions can be placed [3j by 
displaying the NP contribution only, i.e. as if there were no SM contribution at all. This 
procedure is similar to what was traditionally done in the studies of NP contributions to 
mixing, so we shall employ it here too. 

The purpose of this paper is to revisit the problem of the NP contribution to y^ and pro- 
vide constraints on R-parity- violating supersymmetric (SUSY) models as a primary example. 
It has been recently argued in 23J that within R- SUSY models, new physics contribution to 
y-£> is rather small, mainly because of stringent constraints on the relevant pair products of 
RPV coupling constants. However, this result has been derived neglecting the transforma- 
tion of these couplings from the weak isospin basis to the quark mass basis. This approach 
seems to be quite reasonable for the scenarios with the baryonic number violation. However, 
in the scenarios with the leptonic number violation, transformation of the RPV couplings 
from the weak eigenbasis to the quark mass eigenbasis turns to be crucial, when applying 
the existing phenomenological constraints on these couplings. 

We show in the present paper that within R-parity-breaking supersymmetric models with 
the leptonic number violation, new physics contribution to the lifetime difference in D° — D 
mixing may be large, due to the last term in eq. (11.71) . When being large, it is negative (if 
neglecting CP- violation) , i.e. opposite in sign to what is implied by the recent experimental 
evidence for D° — D° mixing. 

The paper is organized as follows. In Section 2 we discuss the R-parity violating interac- 
tions that, in particular, contribute to D° — D° lifetime difference. We confront the form of 
these interactions in the weak isospin basis to that in the quark mass basis, emphasizing the 
important differences. In Section 3 we re-derive formulae for the RPV SUSY contribution 



to y-Q. Unlike ref. 



231 ]. transformation of the RPV coupling constants from the weak to 



the quark mass eigenbasis is taken into account. Also the behavior of different R- SUSY 
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contributions in the limit of the flavor SU(3) symmetry is discussed in details. In Section 4 
we examine the existing phenomenological constraints on the RPV coupling constants. The 
importance of taking into account the transformation of these couplings from the weak to 
the mass eigenbasis is emphasized again. We present our numerical results in Section 5. We 
conclude in Section 6. Appendices contain some details of derivation of bounds on the pair 
products of RPV couplings, relevant for our analysis. 



II. R-PARITY BREAKING INTERACTIONS: WEAK VS MASS EIGENBASES 



We consider a general low-energy supersymmetric scenario with no assumptions made on 
a SUSY breaking mechanism at the unification scales (~ (10 16 — lO 18 )^^). The most 
general Yukawa superpotential for an explicitly broken R-parity supersymmetric theory is 
given by 

(2.1) 



W $ = H ^iikULjEl + X' ijk LiQjDl + -X( jk U?DjD c k 



where Li, Qj are SU(2)l weak isodoublet lepton and quark superfields, respectively; Ef, Uf, 
D\ are SU(2) singlet charged lepton, up- and down-quark superfields, respectively; X^ and 
X[j k are lepton number violating Yukawa couplings, and A" jfc is a baryon number violating 

—Xjik, X'[ jk = —X'La- To avoid rapid proton decay, we assume that 



Hjk ~ an d work with a lepton number violating R- SUSY model. 



Yukawa coupling; 

For meson-to-antimeson oscillation processes, to the lowest order in the perturbation 
theory, only the second term of ( 12. ID is of the importance. The relevant R-parity breaking 
part of the Lagrangian is the following: 



\jk 



—p- d w m w —Ti w d w p- — d w ~ c 



kit. iR JL 



+ h.c. 



(2.2) 



The superscript w indicates that the quark and squark states in (12. 2ft are weak isospin 
eigenstates. The weak and mass quark eigenstates are related by the unitary transformations 
(we assume that left- and right-chiral quarks have the same transformation matrices): 



Uj — S u - n u r 



Sd km d m 



(2.3) 
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where 



and 



^Su^^YujtjSujn — Snn'hun, S*d k , m Xd k , k Sd km ~ ^mm'hd m (2.4) 

3d' k i k ' 



T.S: k] S dkn = V 3n (2.5) 

k 

In ( 12. 4p Y u , Yd are quark-Higgs-quark R-parity conserving Yukawa couplings in the weak 
isospin basis and h u , h d are these couplings in the quark mass eigenbasis. In (12.51) . Vj n 
stands as usually for the (Standard Model) CKM matrix. 

Generally speaking, squark transformation matrices from the weak to the mass eigenstates 
are different from those for quarks. Nevertheless, we choose for squarks to be rotated by the 
same matrices S u and S d that make quark mass matrices diagonal, i.e. 

t.w _ a ~ ~w _ a ~ 

3L Lj u jn u 'n L , u j R tJ v,j n u nR 

d™ L = Sd km d mL , d% R = Sd km d mR (2.6) 

This is a super-CKM basis, in which the squark mass matrices are non-diagonal and result 
in mass insertions that change the squark flavors ^-j^. This source of flavor violation is 
very important in the pure MSSM sector. In particular, it plays crucial role in examining 
the MSSM contribution to D° - D° mass difference 

In the R-parity breaking part of SUSY Lagrangian, flavor changing neutral currents are 
present a priori. In order to simplify our analysis, we put all the squark masses to be nearly 
equal. Then the squark mass matrix is proportional to the identity matrix, i.e. it is diagonal 
in any basis. 

Using (}2~3]) . (12~5|) and (F2J|, one may rewrite as 



&i L d mR u rL + u TL d mR ei L + d mR e iR u rL 



+ 



i,j,k,m,n,r 

+ X! KjkSd km Sdj„ Vi L d mR d nL + d nL d mR Vi L + d* mR v c iR d nL + h.c. (2.7) 

i ,j,k,m,n 

At this point one may redefine, without loss of generality, the couplings A' as 

^ijk^d km ^d Jn — ► \nm (2-8) 

This is also equivalent to choosing the weak and mass eigenbases for down-quarks being the 
same, while for up-quarks they are related by the CKM matrix 2 . 



This definition of A' is not unique. For example, Allanach et al. [28| used the up-quark weak and mass 
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Defining A • 



V* n X' inm and renaming the summation indices, we rewrite (12.71) as 



A* = _ E ^ijfc [ e *L^j? M iL + u j L d k R e iL + d* kR e 

i,3,k 

+ E Kjk [viJk R djL + d jL d kR Vi L + d* kR v c iR d jL 
i,j,k 



"i R U jL 



+ 



+ h.c. 



(2.9) 



As it follows from (12.91) . (s)down-down-(s)neutrino vertices have the weak eigenbasis cou- 
plings A', while charged (s)lepton-(s)down-(s)up vertices have the up quark mass eigenbasis 
couplings A'. 

Very often in the literature (see e.g. {4], j^, [30}- [3^]) one neglects the difference between 
A' and A', based on the fact that diagonal elements of the CKM matrix dominate over non- 
diagonal ones, i.e. 

V jn = 6 jn + 0(\) so \ ijk n\' ijk + 0(\) (2.10) 

where A = sin# c ~ 0.2, with 9 C being the Cabibbo angle. 

Notice that relation Eq. (12.101) is valid if only there is no hierarchy in couplings A'. On 
the other hand, the existing strong bounds on pair products A' x A' (or A' x A') 28|, |30|, [31] 



and relatively loose bounds on individual couplings A' [28] suggest that such a hierarchy may 
exist. As we will see in Section 4, pair products A' x A' may be orders of magnitude greater 
than corresponding products A' x A'. 

To the end of this section, we explicitly write down the terms of the R-parity breaking 
part of the Lagrangian that contribute to D° — D° lifetime difference: 



C 1 



E 



j/ 1 -75\ r, ~ -/!-75 



~* - /I + 75 



c + 



+A,:n e,- it 



d + \%e* r u 



~* _ /I + 75 



E 

k 



\\oi.dl„e c ( — ^] c + X' 22k d* k // f 1 75 ) < 



,/* , - f l + 75 A _ c , w* j _- f l + 75 



+Kikdk R u 



eC + X 21k d ka u 



(2.11) 



2 J Z1K nK \ 2 

In the next section we will integrate out heavy degrees of freedom in (I2.1ip . thus finding 



eigenbases to be the same, relating the bases for down-quarks by the CKM matrix. Another possibility 
is to redefine A' in such a way that (s)up-(s)down-charged (s)lepton vertices have the couplings A' while 
(s)down-down-(s)neutrino vertices have the couplings A' • Vckm 29(. Clearly all these approaches are 
equivalent. 
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d, s 



I 




d, s 
a) 



FIG. 1: D° — D° mixing diagrams with R-parity breaking interactions: a) within the full electroweak theory; 
b) within the low-energy effective theory. In these diagrams, D° - D° oscillations occur via two subsequent 
AC = 1 transitions with the exchange of W boson and a charged "left" slepton, i=l,2,3. 

^,-SUSY part of AC = 1 effective Hamiltonian. Then we will compute R-parity breaking 
SUSY contribution to ATp- 

III. D° - D° LIFETIME DIFFERENCE 

Assuming CP-conservation, the normalized D° — D° lifetime difference is given by 
f 



Ud = 1™ 

2m D T D 



(3.1) 



where H^ c=1 is an effective Hamiltonian including both SM and NP parts. To the lowest 
order in the perturbation theory, R-SUSY contribution to D° — D° mixing comes from the 
one-loop graphs with 

• W ± boson, charged slepton and two down-type quarks (Fig. dh); 

• two charged sleptons and two down- type quarks (Fig. [2k.); 

• two down-type squarks and two charged leptons 3 (Fig. [3^) . 

Within the low-energy effective theory, D° — D° lifetime difference occurs as a result of a 
bi-local transition with two AC = 1 effective vertices. The relevant low-energy diagrams 
in Fig.'s lb) - 3b) are derived by integrating out of heavy boson, charged slepton and 
down-type squark degrees of freedom. 



3 As it follows from (|2.1ip . lepton propagators in Fig. [3] must be constructed by contractions of charge 
conjugates of the electron and/or muon field operators. 



s 



d, s 



c/, s 



a; 




FIG. 2: Same as in Fig. [IJ however both of AC = 1 transitions are due to a charged slepton 
exchange now, i = 1, 2, 3, i' = 1, 2, 3. Both of the effective AC = 1 vertices are NP vertices. 



For R-parity-violating SUSY models one can therefore write 



AC=1 



rrAC=l j_ ttAC=1 , ttAC=1 

n w SM + n w~, + -"vy„- 



(3.2) 



The first term in the r.h.s of (13.21) is the Standard Model contribution, whereas the second 
term comes from AC = 1 transitions with a slepton exchange and the last term comes from 
AC = 1 transitions with a squark exchange. The Standard model part of AC = 1 effective 
Hamiltonian is given by 



ttAC=1 

n w SM 



Gi 



CiOO 5 aiai 5 asa2 + C 2 (/i c ) 5 aia2 5 a ' iai 



V2 

x E V^V^ u a *(x)Y(l - 7 5 )gr(^) ?2 3 (^)7m(1 - 7 5 )c a4 (x) 



(3.3) 



91. 92 



where q\ = s, d, g 2 = s, d, aj are the color indices, and C\ and C 2 are the operator Wilson 
coefficients. The Wilson coefficients are to be evaluated at a low-energy scale /i c , which we 
choose here as \i c = m c . 

To simplify the following calculations, let us assume that all the sleptons and all squarks 
are nearly degenerate, i.e. 



m~ e% = = and m d ~ = = rriq. 



(3.4) 



With this assumption, the low energy effective Hamiltonian for the R-parity-violating inter- 
actions are given by 



H 



AC=1 



X 



Ci(/x c ) 5 aiai 6 asa2 + C 2 (/x c ) $ aia2 5 aaaA 
J2 ^if u<*(x)(l+ l5 )q?(x) g-3(x)(l- 75 )c^(x) 

91, 92 4m £ 



(3.5) 



a) 



b) 



FIG. 3: Same as in Fig.'s [H [2j however both of AC = 1 transitions occur due to exchange of 
down-type squarks now, k = 1,2,3, k' = 1,2,3. Subsequently the intermediate charmless states 
are charged (anti)lepton states. 



and 



H, 



AC=1 



A, 



E t4 *"(*)(! + 75») ^(x)(l - 75 )c a (x) 



(3.6) 



where q\ = s, d, = s, d, t\ = e, /x, and £2 = e, /i. The superscript "c" stands for charge 
conjugation. Also, 



(3.7) 



We assume that A gig2 and A^ 2 are real. 

The insertions of Hamiltonians of eqs. (13.31) . (13. 5p . and (13.61) can lead to the lifetime 
difference in D° — D° system. Let us write it as 



vd = vsm + vsm,np + Vei + Vmi 



(3.8) 



where 



USM,NP 



Im 



D 



mld'xT^H^{x)H^=^ + 



(3.9) 



is the term coming form the interference of the SM and NP contributions to H^F 1 , and 

1 



2m D T D 
1 



Im 



Im 



D 



(D°\i J d 4 x t\h^ c =\x)H^=\^\D q ) , (3.10) 

(3.11) 



(Z) |z / d 4 x T< H^ =1 (x)H^ =1 (0) >\D°) 
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are coming from two insertions of the NP vertices. 

It might seem unreasonable to include double insertions of the NP Hamiltonian to com- 
pute yo, as each insertion generates a contribution that is suppressed by some NP scale 
Mtvp, which in general is greater than the electroweak scale set here by My/. Yet, as the 
Standard Model contribution is zero in the flavor SU(3) limit (i.e. suppressed by powers of 



strange quark mass), New Physics contributions can be large [4J]. Also, as can be seen from 



refs. [4] and 23| . usm,np resulting from the single insertion of the NP Hamiltonian is for- 
bidden in the SU(3) flavor symmetry limit. Thus, double insertion of the NP Hamiltonian 
can be important, especially if this contribution does not vanish in the SU(3) limit! This 
construction can give numerically large contribution to yr> if (M w /M NP ) 2 > (m s /m c ) 2 . 

Note that contribution to Ar^ is nonzero if the intermediate states are the on-mass-shell 
real physical states. It is therefore easy to see from the energy-momentum conservation 
that diagrams like those in Fig.'s [IH but with b-quarks, rr, r/i pairs running a loop, are 
irrelevant for our analysis. While the diagrams with a re pair running in a loop do give 
nonzero contribution to Ar^, their contributions are suppressed by the available phase 
space. Thus, we shall not consider them too. 

It is known that correlation function in (13. ip (as well as those in (I3.9p - fl3.lip ) may be 
presented as a sum of local AC = 2 operators, which corresponds to 1/ m c power expansion 
of (13. ip (or (13.91) - (13. lip ). Here we are interested in the lowest order terms in this expansion. 
Keeping only the leading terms in x s = m 2 s /m 2 c and Xd = m 2 d /m 2 c , we get 



USM,NP 



y/2 A-Km D T D \mjj 
+ A (X ss x s - X dd x d ) - X 2 Xds\ / x s x d 



X s dV X s X d + 

(Q) 



(3.12) 



and 



m 2 c (X 2 ss + X 2 dd + 2X sd Xd s ) 



K 2 



(Q) + 



+ 



K 2 -K x (Q s ) 



where A = sin 6c is the Wolfenstein parameter, and 

■1-75 



(Q) = (D u | n ai (0)7 



c ai (0) w a2 (0hJ^-^) C a2 (0) I/) ) 



(Qs) = (D°\ « ai (0) c ai (0) « a2 (0) 



" ' 2 

'1+75 



c a2 (0) \D°) 



(3.13) 

(3.14) 
(3.15) 



11 



are the matrix elements of the effective low energy AC = 2 operators and 

K x = 3 d C x + d C 2 + C 2 C h K 2 = C 2 C 2 (3.16) 
K x = 3 Of + 2 C x C 2 , K 2 = C\ (3.17) 

are the Wilson coefficients. It is important to stress that ]Jsm,np, just like a Standard Model 
contribution, vanishes in the limit of exact flavor 577(3) symmetry - it is proportional to 
light quark masses via x s , Xd and ^/x s Xd- On the contrary, is nonzero even in the limit 
of exact flavor SU(3) symmetry! Therefore, as we shall see in Section 5, y^ dominates over 
ysM,NP if R-parity breaking coupling products \ ss and/or A^ approach their boundaries. In 
other words, contribution of diagrams in Fig. 2 with both of AC = 1 vertices generated by 
new physics interactions, dominates over the contribution of diagrams in Fig. 1, with one of 
the AC = 1 vertices coming from the Standard Model and the other one coming from new 
physics. 

Similarly, keeping only the leading order terms in x e = m 2 /m 2 , x M = m 2 ^/m 2 c , one gets 



As one can see from (I3.18p . yqq is non- vanishing in the limit of exact flavor SU(3) symmetry 
as well. 

As usual, we parameterize matrix elements (Q) and (Q s ) in terms of B-factors {3J], i.e. 

(Q) = \f D ml B Dt (Qs) = ~ fl ml B S D (3.19) 



where 



B S D = B S D (3.20) 



We shall follow the approach of ref. [l] and neglect QCD running of the local AC = 1 
operators generated by NP interactions. Thus, G\ — and G 2 — 1, or 

K 1 = C 1 (m c ), K 2 = C 2 (m c ), K x = 0, K 2 = 1. (3.21) 

Using (13.191) and ( 13.21 j) . one may rewrite (13.1 2D . (13.131) and ( 13.18P in a following form: 

ysM,NP = — 7=r ^ D ^ mD f — | ) [Ci(m c ) + C 2 (m c )] \\ sd y/x s x d + 
V2 67rr D \mjj L 

+ A (\ ss x s - XddXd) - \ 2 XdsVx^c~d (3.22) 



m 2 c 
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vn 



vm 



- m 2 c f 2 D B D m D 
288vrr D m| 

_ m 2 c fj)B D m D 
2887rr D ml 



2 8 5 D 



hB s D 



8B 



- 1 



(3.23) 



(3.24) 



Formulae fl3.22p - fl3.24p involve only the lowest order short- distance (perturbative) contri- 
bution to D° — D° lifetime difference. Yet, it has been mentioned already that long-distance 
effects play very important role in D° — D° oscillations. In particular, in the Standard Model, 
where the short- distance contribution to has a suppressing factor ~ m s jm c ll|], long dis- 
tance contribution to D° — D° lifetime difference dominates [10]. However, within /£-SUSY 
models we have a different situation. As it is mentioned above, new physics contribution to 
Ud is non- vanishing in the exact flavor SU(3) limit, thus there is no suppression in powers 
of m s /m c in the dominant short-distance NP terms. In what follows, long distance effects, 
which may be interpreted as K^cdI^c power corrections, are subdominant. Thus, they may 
be neglected to the leading-order approximation that is used throughout our paper. 

Further analysis depends on bounds on R-parity breaking coupling constants, so in the 
next section we discuss the existing constraints on these couplings. 



IV. PRESENT BOUNDS ON R-PARITY BREAKING COUPLING CONSTANTS 



Bounds on the R-parity violating couplings A' have been widely discussed in the literature 



28| - [451 ] . Summary of bounds on A^- fc may be found e.g. in [23]. More recent (updated) 
bounds on some A' x A' pair products, coming from the studies of K° — K° and B° — B° 
7i + uu decays, are presented in 



30 



32] and [33j respectively. 



mixing and K + 

It is interesting to note that bounds on RP V couplings coming from K° — K° and B° — B° 
mixing and empirical individual bounds on couplings A^- fc are derived neglecting the difference 
between A' and A'. While for the individual bounds it is a self-consistent approach, for the 
constraints on RPV coupling pair products such an approach in general is questionable. 

Empirical individual bounds on RPV couplings are derived, assuming that only one 
coupling X' i j k is nonzero at a time. If such an assumption is made, then it is easy to see that 



K jk = Kik>< (l + 0(A 2 = sin 2 #, 



c, 



(4.1) 



Knk — 0(\) X X ijk 



(4.2) 
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if n 7^ j, and 

Knm — (4-3) 

if r ^ i or m ^ k. 

Thus, as it follows from fl4.ip - fl4.3l) . when deriving an individual bound on A^ fc by studying 
a given process, there is no essential difference whether the R-SUSY diagram for this process 
contains X' i j k or it contains A- jfc at the vertices. 

Of course, in the realistic R-SUSY scenarios several A' couplings are in general non-zero. 



As it has been pointed out in 28|, even if at the unification scales (~ (10 16 — 10 18 )GeV) one 
has only one non-zero RPV coupling, other non-zero RPV couplings appear when evolving 
down from the unification scales to the electroweak breaking scale. However, the individual 
bounds on A' couplings are still approximately valid, if one assumes that one RPV coupling 
dominates over all other ones. If several couplings dominate, individual bounds may still be 
used, if they are not correlated or weakly correlated with each other. 

The situation with the constraints on the RPV coupling pair products is more compli- 
cated. As we will see, bounds on A' x A' and the corresponding A' x A' products may be 
different by several orders of magnitude. One must therefore be careful when using the 
bounds given in the literature and specify whether these bounds are on A' x A' product or 
they are on A' x A'. This may be easily done, using the following "rule of thumb": 

• If the process that is used to put constraints on the RPV coupling products is described 
by diagram(s) with down-down-sneutrino or down-sdown-neutrino vertices, bounds are 
derived on a A' x A' product. 

• If such a process is described by diagram(s) with up-down-charged slepton, up-sdown- 
charged lepton or sup-down-charged lepton vertices, bounds are derived on a A' x A' 
product. 

• If both types of vertices are present, bounds are derived on some admixture of A' x A' 
and A' x A' products. 

In addition to the individual bounds, we use here constraints on the RPV coupling pair 
products that are derived from study of K + — > -n^vv decay and Am^-o. R-parity breaking 
SUSY contribution to K + — > 7t + uj? is described by tree-level diagrams with a down-type 



squark exchange and quark-squark- neutrino interaction vertices [29|, |33j, |34j. Thus, this 
decay gives bounds on A' x A' products. 
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The situation with K° — K° mixing is more involved: there are several sets of R-SUSY 
diagrams that contribute to this process. In order to get bounds on the RPV couplings, one 
assumes that only a given RPV coupling product or a given sum of RPV coupling products 
is nonzero. Possible bounds on the RPV coupling pair products have been originally listed 
in 31] . Recently these bounds have been improved in [30J . Bounds that are relevant for our 
analysis are presented in Appendix A. We also specify which of them are for A' x A' pair 
products and which of them are for A' x A'. 

Keeping in mind everything that has been said above, let us consider the RPV coupling 
products, which are present in formulae (I3.22p -( l3~24l) . We start with 

A ss = ^2 K*uK22 = VlnVty X'* n2 X' ij2 . (4-4) 

i i,j,n 

Using Wolfenstein parametrization for the CKM matrix, keeping for each A' x A'* product 
only the leading order term in A = sin 8c, and assuming that all A' x A'* products are real 
(no new source of CP- violation), we rewrite (14. 4p in a following form: 

Ass = Y K*i2 K22 — Y K*u K22 + ^[H l-^l 2 - F IKul 2 

i i i i 

+ AX 2 £ Af 12 A:3 2 + AX 3 (1 + P -i V )j: X% 2 K 22 

i i 

+A 2 A 5 (p-^)£|Ay 2 (4.5) 

i 

There is a strong bound on the Cabibbo-favored term in the r.h.s. of (14. 5p from the K + — > 



-n + vi> decay. Assuming that X'* lk X' i2k 7^ only for k=2, one gets [33j 



i A - A -i S6 - 3 - 10 ' 5 (3o^) <«» 



We have rescaled the bound of ref. 



33| | to the units of m 9 -/ '300 GeV. Values of the squark 



masses less than 300 GeV are disfavored by many experiments (see 46[ for more details). 



For this reason, we follow ref. 30] assuming that rrig > 300 GeV. 

If squarks happen to be superheavy 4 , there is still a strong bound on the Cabibbo favored 
term in ( 14.51) from K° — K° mixing. As it follows from our discussion in Appendix A, 

IF A; 22 | < 2.7 x 10- 3 ( m J V (4.7) 
112 % 22\ - VlOOGeW v ; 



We thank X. Tata for discussion of this scenario. 
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Thus, the Cabibbo favored term in (I4.5p is strongly suppressed, if one assumes that only 
K12 0- ari d K22 0- On the other hand, even under such an assumption, one still has 



Ass = K12 K22 7^ K12 K 
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due to the first order Cabibbo suppressed terms in (14.51) . Furthermore, constraints (14. 6 p 
or (14. 7\\ may in particular be satisfied, when |A^ 22 | is close to its boundary value whereas 
|A- 12 | — > 0, and vice versa. Taking into account that individual bounds are, in general, orders 
of magnitude looser than (14.61) or (14. 7p . it is not hard to see that X ss is dominated by the 
first order Cabibbo suppressed term in (14. 5p . 

Further on we will very often deal with a situation, when expanding A' x A' products in a 
basis of A' couplings, the Cabibbo favored term is negligible whereas the first order Cabibbo 
suppressed term dominates, and the only possible constraints on the first order Cabibbo 
suppressed term are the individual bounds on A' couplings. In order to use these bounds we 
assume hereafter that only one coupling A^- fc dominates at a time. 

After making such an assumption, it is easy to see that 

30Q J eV j < Ass < 0.29, if m- q <\TeV, 

-0.29 < X ss < 0.29, if m- q > ITeV (4.8) 

The upper bound on X ss is derived when one of A^ 22 couplings dominates. Individual bounds 



on X' i22 are the loosest for i = 3 [28| . For rrig > 300GeV, | A322 1 < 1.12 - this is the pertur- 
bativity bound on A322- The lower bound on X ss is derived when one of X' il2 couplings dom- 
inates. Individual bounds on X' il2 are the loosest for i=3 again: |A' 312 | < 0.33(mq/300GeV^), 
if rrig < ITeV and |A 3 i 2 | < 1.12 - the perturbativity bound, if rrig > ITeV. 
It is important to stress that, in general, as it follows from ( I4.6p . (I4.7p . ( 14 .81) . 

A ss = E^ 2 4 »EAf 12 ^ 22 (4.9) 

i i 

Thus, as it has been already pointed out in the beginning of this section, bounds on A' x A' 
products differ by several orders of magnitude from those on corresponding A' x A' products. 
In the considered case, A' x A' product is restricted by much weaker bound than corresponding 
A' x A' product. 

Relation (14.91) plays crucial role in our analysis. We will see in the next section that, as 
a consequence of this relation, R-parity breaking SUSY contribution to Ar^ is quite large. 
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For Xddj analysis is performed in exactly the same way and yields 

- 0.025 -f— < X dd < 0.29, if ma < ITeV, 



\300GeW 

-0.29 < X dd < 0.29, if rriq > ITeV (4.10) 
Also, the relation similar to (14. 9p is obtained: 

\« = £\^>E*Si*m (4.11) 

j i 

and relation (14.111) is as crucial as (14.91) . It is also useful to transform (14 .8p and (14.101) onto 
restrictions on A 2 S and X dd : 

AL « A 2 [£ |Ay 2 - £ |A; 12 | 2 ] 2 < 0.0841 (4.12) 

i i 

AL « A 2 [E |A; 21 | 2 - £ |A: n | 2 ] 2 < 0.0841 (4.13) 

i i 

Bounds on X ds and X sd are derived using the experimental data for Am^o. As it follows 
from formula (1A.ll) in Appendix A, 

M-|E^A; 22 |<1.7.10-( I5 ^ F ) 2 (4.14) 

In order to derive constraints on X sd , one must write it in a following form (using A^ 



ijk 



VnjKnk) 



A S d = Aii 2 Ai2i - (^11^22) 1 



(4.15) 



Y A S2A m - Y VjlVnZ ( Y Kj2Knl / 
i Li j,n i 

where prime indicates that the sum over j and n does not contain the term with j — 1 and 
n = 2. Bounds on the terms present in r.h.s. of (14.151) are given in Appendix A. Using these 
bounds, one can see that 

X sd < few x 10- 7 I m J V (4.16) 
VIOOGeW V ; 

It is interesting to note that such strong constraints on X ds and on X sd are derived assuming 
that only one A' x A' or A' x A' product is nonzero. It is also assumed that pure MSSM 
sector gives negligible contribution to Am K o (3^]. These two assumptions are not necessarily 
true. If one gives up these assumption, then destructive interference of the pure MSSM and 
Fi-SUSY diagrams or the one of different JK.-SUSY diagrams will somehow distort bounds 
(I4.15p . (14.161) . However, unless there is a fine-tuning or an exact cancelation between two 
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(or more) diagram contributions, it is very unlikely for the distortion of these bounds to be 
such that Xd s and/or X sd be ~ 10 _1 or ~ 10~ 2 . Therefore in our numerical calculations we 
will use the following relations: 

X ds <^X ss ,X dd (4.17) 
X sd <^X ss ,X dd (4.18) 

For the remaining four coupling products - A ee , A w , A^ e and X efl - that are contained in the 
expression (3.26) for the analysis is similar to that for X ss and A^. For the details and 
subtleties of the analysis, we refer the reader to Appendix B. Here we only point out that 
bounds on A ee , A MAt are the following: 

2 / 



2 / /vyi _ \ 2 



(171 \ f 771 ~ \ 

mckv) £ A » 5 °' 091 (mkv) • if m * £ 530 GeV - < 419 > 

(Tfl - \ ^ 
— < X uu < 0.29, if m g > 530 GeV. (4.20) 
300GeW ~ w ~ 9 ~ V ; 

Also, for two other couplings we get 

/ Tfl \ 2 / Tfl ~ \ 2 

/ 777 \ / r ji ~ \ 

|A M J< 0.033 -I— , |A e J< 0.033 -|— , if m„- > 530 GeV (4.21) 

1 ~ V300GeW ' 1 Ml ~ V300GeW ' q ~ K ' 



We also obtain that 



A Me « A eM (4.22) 



As increases, squark mass dependent empirical bounds on the RPV couplings are replaced 
by squark mass independent perturbativity bounds. In formulae (I4.19p - fl4.2ip . we indicate 
the change in the behavior of the bounds with the squark mass, if it occurs for rrig < ITeV. 
When transforming ( 14. 19f) -( f4T22l) onto the restrictions on A 2 e , A 2 ^, A^A^, one can see 



that these restrictions are much weaker than the relevant constraints listed in ref. 23J]. This 
is because in the present paper we do not neglect the transformations of RPV couplings 
from the weak eigenbasis to the quark mass eigenbasis. More precisely, we do not neglect 
the difference between A' x A' and A' x A' pair products. 

From fl4.19p - fl4.22p . one can also see that generally speaking, 

AJ M » A /te A eM » A 2 e (4.23) 
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It is worth mentioning here that additional bounds on A ee , A w , A Me , X efJi may be derived 
from studying rare D- meson decays, such as D — » X£ + £~, D° — > £ + £~, etc [2l| . As it 
follows from the analysis performed in ref. 21], bounds derived in this way may be even 
stronger than those given by (14.191) - (|4.2l|) . Bounds coming from the rare D- meson decays are 
however still to be elaborated in details, taking into account new experimental data, as well 
as possible impact of the long-distance SM and (short-distance) pure MSSM contributions. 
Such an elaboration is beyond the scope of this paper, in particular because turns to 
be a (numerically) subdominant part of the new physics contribution to D° — D° lifetime 
difference, even if we use constraints on A ee , A MM , A^e, X efl given by fl4.19n - fl4.2in (see the 
next section). 

Having obtained constraints on all RPV coupling products in (I3.22p - (13.241) . we may pro- 
ceed to computation of Usm,np, Va-, Vqq- 



V. NUMERICAL ANALYSIS 



In our numerical calculations we use [46| Gf = 1.166 ■ 10 5 GeV 2 , A ~ 0.23, Tn 
1.6 • 10~ 12 GeV, m D « 1.865 GeV; m c = m c (m c ) « 1.25 GeV, m s (2GeV) « 95 MeV, 

12/25 21 \ 

m s [m c ) 



m s (m c ) ~ m s (2GeV) 



l a s {m c ) 
\a s (2GeV)) 

1.208 M, B D w 0.; 



105 MeV, x s 



ml(m c ) 



0.007; 



11 



D 



0.22 



C x (m c ) = -0.411, C 2 (m c ] 

While the value of Bp is known from the lattice QCD calculations, there is no theoretical 
or experimental prediction on Here we follow the approach of ref. ll| . assuming that 



B b D = B 



Di B D = 0.8B D , 



Bjj = 1.2B D . 



(5.1) 



Let us first determine the sign of i/sm,np, Vlh Vqq- Using relations (I4.17p . (I4.18p . (I4.23p . one 
may rewrite equations (I3.22I) - (I3.24I) in a much simpler form, 

- G F flB D m D ( m 2 c 



USM,NP 



Vie 



V2 QnT D 
- m 2 c fj)B D m D 
288vrr D m\ 

ml j 2 D B D m D 
2887rr D mj 



m 



CAm c ) + C 2 (m c 



A Ace Xc 



5B 
8B 



2 8B D 



- 1 



D 



It follows from (I5.2p . (I5.3P that the sign of Usm,np is opposite to that of X ss and 



(5.2) 
(5.3) 
(5.4) 
< 0. 
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One can see from (15.41) that the sign of yqq is determined by the factor 



. As 



it follows from (13.201) and (15.11) . for m c = m c (m c ) ~ 1.25GeV, this factor is positive, hence 

Vqq > 0. 

and hence yqq flips its sign when using the charm quark pole 



On the other hand 



B D 1 



mass 5 , mP c o[e « 1.65 GeV. 

In general, such an ambiguity in sign of yqq may cause a trouble in numerical evaluation of 
the results, signaling the need for next-to-leading order evaluation of the appropriate contri- 
butions, where the scheme ambiguity cancels out. Here we disregard this sign ambiguity, as 
yqq turns to be a (numerically) subdominant part of the new physics contribution to D° — D° 



lifetime difference. In our opinion, the use of the MS charm mass, m c (m c ) = 1.25 GeV, is 
more appropriate in this calculation. Then yqq has positive sign. 

Let us proceed to our results. It is convenient to start with yqq. Using the listed numerical 
values of parameters present in (15.41) . we get 

/300GW\ 4 



B D = 0.8B D : yqq « 0.0011 A 



/'/' 



B S D = B D : yqq « 0.0038 A> (2°^V (5.5) 



B S D = 1.2B D : yqq « 0.0064 \\ 



mi 



V m 9 / 

/300GWy 



As it follows from (15. 5p . to the lowest order in the perturbation theory, yqq is highly sensitive 



to the choice of parameters B^ and Bd- Moreover, if one uses the approach of ref. 23], 
choosing = Bd or B D = (m 2 c /m D )BD ~ 0.45-Bd, yqq hips the sign 6 . 
Using the bounds on A w given by (14.201) yields 



B S D = 0.8B D : yqq < 0.9 ■ lO" 5 
B S D = B D : yqq < 3.12 -10- 5 (5.6) 
B S D = 1.2B D : yqq < 5.34 ■ lO" 5 



5 To derive the proper value of vnP° le , the two-loop relation between the pole and M S quark masses must 
be used. This is because the MS value of the c-quark mass has been extracted using the perturbative 
QCD analysis up to the order a 2 s [46| . One can check that the use of the three loop relation between the 
pole and MS quark masses [4^] leads to the physically meaningless result m v ° le w 1.93 GeV > mo. 

6 y^q is equivalent to —y(Rpv-RPV,l) m the notations of [23j. 
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for nig < 530 GeV and 



B S D = 0.8B D : ^<0.9.10- 5 (^^y 



B b D = B D : ^<3.12-10' 5 (5.7) 

V m « / 

B S D = 1.2B D : y« < 5.34 • lO" 5 (^^Y 

for nig > 530 GeV. 

Thus, if using bounds on A ee , A MM , A^e, A e/Lt , given by (14.191) - (I4.22p . one obtains that 
Uqq is at least by two orders of magnitude less than the experimental value of ?/d- As it 
was mentioned above, constraints on A ee , A MM , A Me , A eA1 and hence on may become even 
stronger if one elaborates the constraints on RPV couplings coming from the rare D-meson 
decays. Further on we simply disrespect x/qq because of its smallness. This way we also 
avoid the problems related to the dependence of the obtained results on the choice of the 
renormalization scheme and .Bp-factors. 

Consider Usm,np now. For this quantity one gets 

2 

V m i ) 
which after using (14. 8 p yields 



y S M,NP ~ 0.0040 \ ss ( 100GeV \ (5. 8 ) 



for rriq < 1 TeV and 



flOOGeVV flOOGeVV 



0.0011 < y SM ,NP < 0.0011 (5.10) 

V m i J \ m l J 



for rriq > 1 TeV. 



As it follows from (15.91) . (15.101) . \usm,np\ ma Y be by an order of magnitude greater than 



it was quoted in 



23j 7 . This is because the analysis in ref. |23] has been restricted by 



consideration of rriq = 100 GeV only. On the other hand, as it follows from Table I of 



ref. 



28] and our analysis in Section 4, bounds on RPV couplings and hence on \ ss become 



weaker for the greater values of squark masses. Else, unlike ref.'s |4j, [23J], we obtain that 



7 Vsm,np = -V(sm~rpv) in the notations of 23] 
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Vsm,np can be both positive and negative. This is because, as one can see from equation 
(14. 5 p and the following it discussion, X ss may have both of signs even if one assumes that all 
RPV couplings are real and positive. 

Finally, consider y~f v Using the numerical values of the parameters present in f)5.3p . one 

gets 



B b D = 0.8B D : y E « -1.25 



Ai + A 



da 



(lWGeV\ 



4 



B b D = B D : y u » -1.47 X 2 SS + X 2 dd (5.11) 



B b D = 1.2B D : y u « -1.69 



AL + AL 




As one can see from (15.111) . varying the ratio Bp/Bp from 0.8 to 1.2, one gets about 15% 
uncertainty in the predictions for y-f t Thus, yu is only weakly sensitive to the choice of the 
parameter As we are interested in the order of the effect only, we may for a simplicity 
assume £>£ = B D hereafter. 

To be consistent with a one dominant coupling approximation, we will assume that only 
one of the coupling products \ ss or is at its boundary at a time. Notice however that 
if we allow both X ss and Xdd to be simultaneously large, our results will change at most by 
a factor two, which is inessential, if one is interested in the order-of-magnitude of the effect 
only. 

Using the bounds on X 2 S and X 2 dd given by (14.121) and (14. 13j) we obtain 

-0,2p^) 4 , fe -<0 (, 12) 
It is important to stress that \y^\ may be ~ 10 _1 , if = 100 GeV. 

-ii 



This result is in contradiction with the one of ref. |23j: yRpv-pnv,q — ~Uu ^ 2-5 • 10 



for mi = lOOGeV. This contradiction is related to the fact that authors of ref. [23] , following 
other papers on the meson-antimeson mixing phenomenon, have neglected the transforma- 
tion of the RPV couplings from the weak eigenbasis to the quark mass eigenbasis. This 
allowed them to impose very stringent constraints on X 2 S and X dd from K + — > -n^vv decay. 
As it follows from our discussion in Section 4, this approach is not always appropriate 8 . 



Unless one imposes the conditions A- 22 ~ A- 12 and A- 21 ~ A- n . 
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We are now able to compute the total New Physics contribution to D° — D° lifetime 
difference, 

Vnew = VSM,NP + Vil + Vqq- 

As it is mentioned above, we neglect y^ because of its smallness. Also, as it follows from 
(15. 8p and (15.111) . y^ ^> usm,np unless Add — ► and the ratio \ ss /m~ is small enough. It is 
not very hard to see after doing some algebra that 

flOOGeVV , 

-0.12 < yn + y S M,NP< 2.72 • 10~ 6 (5.13) 

V m i ) 

The (negative) lower bound in (15 . 1 3[) is derived neglecting ysM,NP as compared to y~f t The 
(positive) upper bound in (15 . 1 3[) is derived for Add = and X ss = —0.00136 (mjj/lOOGeV) 2 , 
when ysM,NP = —tya- As it follows from (15.61) and (15.131) . y new is negligible, if positive, and 
may be as large as ~ 10 _1 , if negative. 

Thus, within the R-parity breaking supersymmetric models with the lepton number vi- 
olation, new physics contribution to D° — D° lifetime difference is predominantly negative 
and may exceed in absolute value the experimentally allowed interval. In order to avoid a 
contradiction with the experiment, one must either have a large positive contribution from 
the Standard Model, or place severe restrictions on the values of RPV couplings. As it 

n 

follows from [lOj], ysM may be as large as ~ 1%. In what follows, \y n ew\ must be ~ 1% or 
smaller as well. If \y new \ ~ 1%, one may neglect ysu,NP as compared to y^. Then, imposing 
condition 

- 0.01 < y new w y s (5.14) 

one obtains that either > 185GeV, or if < 185GeV, condition (15.141) implies new 
bounds on X ss and Add: 

2 



|A„| < 0.082 



(7TI ~ \ ^ 

Note that bounds (15.151) and (15. 161) may not be saturated simultaneously. (15.151) is saturated 
if Add = 0. Subsequently, (15.161) is saturated if X ss = 0. For the opposite limiting case, 
Ass = Add, one gets v^2 times stronger restrictions: 

mkv) • l A -l 5 °' 058 {Tmkv) (5 - 17) 
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It is interesting to compare the restrictions on A ss and Add, given by (I5.15l) - (l5.17p . with those 
derived in ^] from study of D° — D° mass difference. Translated to our notations, we may 
rewrite the relevant constraints of ref. [3j in the following form: 

X„ < 0.085^ , X M < 0.085^ (^) (5.18) 

This constraint has been derived assuming that rriq = m^. If rriq 7^ m^, bounds in (I5.18P 



must be divided by the factor §y 1 + rn~/m~, as it follows from formulae (130)- (134) of ref. 
jsj]. Assuming for a simplicity that m~ ^> m~ and inserting x exp = 0.0117 into (15.181) . one 



q 

gets 

^ 0037 imkv) (519) 



Ks < 0.0037 

Thus, bounds of on X ss and Add are about 20 times stronger than our ones. On the other 
hand, constraints of ref. [3] on the RPV coupling products are derived in the limit when the 
pure MSSM contribution to Am^ is negligible. Generally speaking, the MSSM contribution 
to D° — D° mass difference is significant even for the squark masses being about 2GeV. 
In what follows, the destructive interference of the pure MSSM and ]&-SUSY contributions 
may distort bounds (15.191) . making them inessential as compared to fl 5 . 1 5 [) - (1 5 . 1 T[) or even to 

d4SD, duo}. 

Contrary to this, pure MSSM contributes to Ar^ only in the next-to-leading order via 
two-loop dipenguin diagrams. Naturally, this contribution is expected to be small. In what 
follows, unlike those of ref. [jj], our constraints on the RPV coupling products X ss and Add, 
given by (I5.15p - (l5.17p . seem to be insensitive or weakly sensitive to assumptions on the pure 
MSSM sector of the theory. 

Thus, our main result is that within the R-parity breaking supersymmetric theories with 
the leptonic number violation, new physics contribution to Ar^ may be quite large and is 
predominantly negative. 

For simplicity we assumed that all sleptons have nearly the same mass and all squarks 
have nearly the same mass. It is easy to see that taking into account the difference between 
the slepton masses does not affect our main results. There are however subtleties concerning 
to the squark masses. First, recall that our analysis has been performed for trig > 300 GeV. 
While this constraint is quite reasonable for d and s, bottom squark is still allowed experi- 



mentally to be about 100 GeV 46]. On the other hand, we have seen that bounds on i/sm,np 



and yn either grow or are insensitive to the squark masses. As for the bound on it is 
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insensitive on rriq for low values of the squark masses. Thus, no new effect is going to be 
observed, if one takes the squark masses to be about lOOGeV. 

Another point to be made, is that the squark mass matrix is in general non-diagonal 
in the super-CKM basis, if one takes the squark masses to be different. In this case, to 
take properly into account the squark mass insertion effects, one should also give up the 
simplifying assumption that left- and right-chiral quarks (of a same flavor) have a same 
transformation matrix from the weak eigenbasis to the mass eigenbasis. It has been already 
mentioned in Section 2, that no new flavor violation effects are obtained, however this may 
somehow weaken bounds ( 14.19P - (14.211) on A ee , A MM , A Me A eM , when applying arguments 
analogous to those used in Section 4. However, as it was mentioned above, A ee , A w , A Me 
A eA1 are expected to get additional strong constraints from the analysis of the rare D-meson 
decays, so that one may expect for y^ to be in any case restricted by even more stringent 
bound than (15. 5p . In other words, giving up the assumption of nearly equal squark masses 
leads to complication of the analysis without observation of any new effect. If being large, 
RPV SUSY contribution to the lifetime difference in D° — D° mixing still may have only 
negative sign. 

When studying the lifetime difference in D° — D° mixing within the Standard Model 
and beyond, one usually assumes that CP- violating effects are negligible {4, 10, 11, 22, 23]. 
Following this strategy, we have chosen for the RPV coupling products that contribute to 
D°-D° mixing amplitude to be real. The natural question arises if our results may be 
affected by possible complex phases of these coupling products. Clearly, \y new \ still may be 
large, however the complex phases may possibly affect its sign. One may suggest - because 
of no evidence of CP-violation in D° — D° system {5, [f| - that the phases of the relevant 
RPV coupling products are small. In this case, contribution to D° — D° lifetime difference, 
proportional to the imaginary parts of the RPV coupling products, is subdominant and 
cannot affect the sign of y new - if being large in the absolute value, y new is negative . Yet, 
it may happen that RPV coupling products that contribute to D° — D° mixing have large 
phases, and no evidence of CP-violation in D° — D° system is related to the fact that - 
unlike the D° — D° oscillations - R-SUSY contribution to D° meson decays is rather small. 
In that case the formalism, used in our paper, is not valid anymore. More general and 
involved approach should be used, taking into account possible correlations in the values of 
D° — D° mass and lifetime differences as well as possible correlations in the SM, pure MSSM 
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and RPV sector contributions. Thus, to clarify if the RPV couplings complex phases may 
affects the sign of the NP contribution to D° — D° lifetime difference, thorough and detailed 
study of the case, when the relevant phases are large, is needed. 

VI. CONCLUSION 

We computed a possible contribution from R-parity- violating SUSY models to the lifetime 
difference in D° — D° mixing. Even though the D° — D° system is rather unique in that 
the Standard Model predicts vanishing of yo in a symmetry limit, the technique and results 
described here can be applied to other heavy flavored systems, especially those where the 
the Standard Model predictions are very small, such as .B^-system. The contribution from 
RPV SUSY models with the leptonic number violation is found to be negative, i.e. opposite 
in sign to what is implied by recent experimental evidence, and possibly quite large, which 
implies stronger constraints on the size of relevant RPV couplings. 

We discussed currently available constraints on those couplings (especially on the products 
of them), available from kaon mixing and rare kaon decays. We emphasize that the use of 
these data in charm mixing has to be done carefully separating the constraints on RPV 
couplings taken in the mass and weak eigenbases, given the gauge and CKM structure of 
D°-D° mixing amplitudes. 
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APPENDIX A: BOUNDS ON THE RPV COUPLING PAIR PRODUCTS FROM 

Am K o 

R-parity breaking part of SUSY contributes to K° — K° mixing by the tree-level diagram 
with a sneutrino exchange, by the so-called L2 type of box diagrams with W ± boson and a 
charged slepton exchange and by the so-called L4 type of box diagrams with all four vertices 



26 



being new physics generated vertices [30|. Bounds on the RPV coupling products are derived 
assuming that only a given pair product or a given sum of pair products is non-zero. 

Here we list the bounds, derived in [30), that are relevant for our analysis. We consider 
only the case when the pair products are real. We specify which of constraints are for A' x A' 
products and which of them are for A' x A': 
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If one assumes that the RPV coupling products are non-zero only for a given i and a given 
k, one may apply them to each term in the above sums. 

Bounds (IA.1I) - ( 1A.5I) are derived from charged slepton mediated L2 diagrams and (1A.6I) 
is derived from a tree level sneutrino mediated diagram. Naturally these bounds scale with 
the slepton mass squared. Contrary to this, to derive (1A.7|) . both sneutrino mediated and 
squark mediated L4 diagrams are used. Thus, it is not easy to scale this bound. However 
for mi = lOOGeV and = 300GeV, the squark mediated diagrams contribution is about 
10% of that of the slepton mediated ones |30|]. In what follows, (1A.7I) is also approximately 
valid if ^> m~ t Then this bound may be scaled with the slepton mass squared as well. 
Assuming that A^ fc A^ 2fc 7^ only for a given value of k, one gets 
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We do not use bounds of [30j for ij2 x ijl combination products. Using our "rule of thumb" 
one can see that these are bounds on some admixture of A'toA'.-i and A'* 9 A'„i . We use instead 
earlier bounds of ref. [3l| . These bounds are derived using L2 diagrams only, neglecting L4 
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ones. These diagrams vertices contain A' couplings, but not A'. Thus one has 
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APPENDIX B: BOUNDS ON A ee , A w , A e/ ,, A Me 
We may present A ee , A MM , A^e, A eM in a following form: 
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The Cabibbo favored terms in (1B.ll) -( lB~4l) have severe constraints e.g. from study of K + 
7t + vu decay [331 ] : 
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For i = i', bounds are about 30% weaker because of the impact of the SM and pure MSSM 
contributions 

It turns out that because of the stringent bounds on the Cabibbo favored terms, r.h.s. of 
(IB.ip - flB.4l) are dominated by the first order Cabibbo suppressed terms. 

The analysis for A ee and A M/J is very similar to that for X ss and Add- Assuming that one 
of the couplings Ai 2 /c or An^ dominates (say for k=3), one gets 
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In analogous way, assuming that one of the couplings \22ie or X 2 ik dominates, one gets 
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The upper bound in the second line of (1B.8P comes from the perturbativity bound on X' 22k for 
k=2,3 [28(: A 22 fc < 1.12. We indicate the perturbativity bound saturation if only it occurs 
for nig < ITeV. 

The analysis for A Me and A eA1 is more subtle: instead of individual couplings squared in 
absolute value, the first order Cabibbo suppressed terms contain RPV coupling pair products 
now. On our knowledge, there is no bounds on pair products 9 X' 12k X 22k an d ^lik^'iik- Thus, 
we must use individual bounds on these four couplings. As we deal with a pair product, 
we may not anymore assume that only one RPV coupling dominates. We must now allow 
for two RPV couplings to be at their boundaries at a time. There is however one subtlety: 
one may do this, if only there is no correlations between the constraints on A' 22fc and X' 12k or 
between those on X' 21k and A' llfc . 

One can check that constraints on X' 22k and X' l2k are indeed independent of each other and 
constraints on X' llk are independent of the values of XL k . The sources of these constraints 
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On the other hand, one can see from Table I in ref. 28j that 
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Thus, condition (IB. 91) is satisfied to a good extent, when X' uk and X 21k are at their boundaries. 
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28| . to get constraints on the pair products A^A^ and A^A^- Using these 



9 One can meet some bounds in the literature on A' lmfe A2* mfe from study \i — > e-j decay (see 45] and references 
therein). However, using our "rule of thumb", it is easy to see that these are bounds on A' 12fc A2* 2fe , thus 
they may not be used here. 
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constraints and assuming that only one of these pairs is non-zero (dominant) and only for a 
given k (say k=3), one gets 
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In deriving (IB.llj) . one must take into account that products A'^A^ and A'^A^ may be 
both positive and negative. 

Coincidence of bounds on A Me and A eAl is not accidental: the first order Cabibbo suppressed 
terms in equations (IB. 31) and (1B.4I) are complex conjugates of each other. Thus, A Me ~ A*^ 
or because we assume that RPV coupling products relevant for our analysis are real, one 
has 

A^ e « A eM (B.12) 

When deriving (IB.llj) and (IB. 120 . we neglected 0(A 2 ) Cabibbo suppressed terms in the 
expressions for X efl and A^ e . If one assumes that two RPV couplings dominate at a time, one 
should take into account these terms as well. We leave for the reader to verify that 0(A 2 ) 
terms in the expressions for X efl and A^ e have at least several times stronger bounds than 
the first order Cabibbo suppressed terms. 
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